A comparative genomics approach reveals a local genetic signature of Leishmania tropica in Morocco.

In Morocco, cutaneous leishmaniasis (CL) caused by Leishmania (L.) tropica is an important health problem. Despite the high incidence of CL in the country, the genomic heterogeneity of these parasites is still incompletely understood. In this study, we sequenced the genomes of 14 Moroccan isolates of L. tropica collected from confirmed cases of CL to investigate their genomic heterogeneity. Comparative genomics analyses were conducted by applying the recently established Genome Instability Pipeline (GIP), which allowed us to conduct phylogenomic and principal components analyses (PCA), and to assess genomic variations at the levels of the karyotype, gene copy number, single nucleotide polymorphisms (SNPs) and small insertions/deletions (INDELs) variants. Read-depth analyses revealed a mostly disomic karyotype, with the exception of the stable tetrasomy of chromosome 31. In contrast, we identified important gene copy number variations across all isolates, which affect known virulence genes and thus were probably selected in the field. SNP-based cluster analysis of the 14 isolates revealed a core group of 12 strains that formed a tight cluster and shared 45.1 % (87 751) of SNPs, as well as two strains (M3015, Ltr_16) that clustered separately from each other and the core group, suggesting the circulation of genetically highly diverse strains in Morocco. Phylogenetic analysis, which compared our 14 L. tropica isolates against 40 published genomes of L. tropica from a diverse array of locations, confirmed the genetic difference of our Moroccan isolates from all other isolates examined. In conclusion, our results indicate potential regional variations in SNP profiles that may differentiate Moroccan L. tropica from other L. tropica strains circulating in endemic countries in the Middle East. Our report paves the way for future research with a larger number of strains that will allow correlation of diverse phenotypes (resistance to treatments, virulence) and origins (geography, host species, year of isolation) to defined genomic signals such as gene copy number variations or SNP profiles that may represent interesting biomarker candidates.

TAD boundary deletion causes PITX2-related cardiac electrical and structural defects.

While 3D chromatin organization in topologically associating domains (TADs) and loops mediating regulatory element-promoter interactions is crucial for tissue-specific gene regulation, the extent of their involvement in human Mendelian disease is largely unknown. Here, we identify 7 families presenting a new cardiac entity associated with a heterozygous deletion of 2 CTCF binding sites on 4q25, inducing TAD fusion and chromatin conformation remodeling. The CTCF binding sites are located in a gene desert at 1 Mb from the Paired-like homeodomain transcription factor 2 gene (PITX2). By introducing the ortholog of the human deletion in the mouse genome, we recapitulate the patient phenotype and characterize an opposite dysregulation of PITX2 expression in the sinoatrial node (ectopic activation) and ventricle (reduction), respectively. Chromatin conformation assay performed in human induced pluripotent stem cell-derived cardiomyocytes harboring the minimal deletion identified in family#1 reveals a conformation remodeling and fusion of TADs. We conclude that TAD remodeling mediated by deletion of CTCF binding sites causes a new autosomal dominant Mendelian cardiac disorder.

Positive selection in the genomes of two Papua New Guinean populations at distinct altitude levels.

Highlanders and lowlanders of Papua New Guinea have faced distinct environmental stress, such as hypoxia and environment-specific pathogen exposure, respectively. In this study, we explored the top genomics regions and the candidate driver SNPs for selection in these two populations using newly sequenced whole-genomes of 54 highlanders and 74 lowlanders. We identified two candidate SNPs under selection - one in highlanders, associated with red blood cell traits and another in lowlanders, which is associated with white blood cell count - both potentially influencing the heart rate of Papua New Guineans in opposite directions. We also observed four candidate driver SNPs that exhibit linkage disequilibrium with an introgressed haplotype, highlighting the need to explore the possibility of adaptive introgression within these populations. This study reveals that the signatures of positive selection in highlanders and lowlanders of Papua New Guinea align closely with the challenges they face, which are specific to their environments.

Upstream open reading frame-introducing variants in patients with primary familial brain calcification.

More than 50% of patients with primary familial brain calcification (PFBC), a rare neurological disorder, remain genetically unexplained. While some causative genes are yet to be identified, variants in non-coding regions of known genes may represent a source of missed diagnoses. We hypothesized that 5'-Untranslated Region (UTR) variants introducing an AUG codon may initiate mRNA translation and result in a loss of function in some of the PFBC genes. After reannotation of exome sequencing data of 113 unrelated PFBC probands, we identified two upstream AUG-introducing variants in the 5'UTR of PDGFB. One, NM_002608.4:c.-373C>G, segregated with PFBC in the family. It was predicted to create an upstream open reading frame (ORF). The other one, NM_002608.4:c.-318C>T, was found in a simplex case. It was predicted to result in an ORF overlapping the natural ORF with a frameshift. In a GFP reporter assay, both variants were associated with a dramatic decrease in GFP levels, and, after restoring the reading frame with the GFP sequence, the c.-318C>T variant was associated with a strong initiation of translation as measured by western blotting. Overall, we found upstream AUG-introducing variants in the 5'UTR of PDGFB in 2/113 (1.7%) undiagnosed PFBC cases. Such variants thus represent a source of putative pathogenic variants.

De Novo Variants Found in Three Distinct Schizophrenia Populations Hit a Common Core Gene Network Related to Microtubule and Actin Cytoskeleton Gene Ontology Classes.

Schizophrenia (SZ) is a heterogeneous and debilitating psychiatric disorder with a strong genetic component. To elucidate functional networks perturbed in schizophrenia, we analysed a large dataset of whole-genome studies that identified SNVs, CNVs, and a multi-stage schizophrenia genome-wide association study. Our analysis identified three subclusters that are interrelated and with small overlaps: GO:0007017~Microtubule-Based Process, GO:00015629~Actin Cytoskeleton, and GO:0007268~SynapticTransmission. We next analysed three distinct trio cohorts of 75 SZ Algerian, 45 SZ French, and 61 SZ Japanese patients. We performed Illumina HiSeq whole-exome sequencing and identified de novo mutations using a Bayesian approach. We validated 88 de novo mutations by Sanger sequencing: 35 in French, 21 in Algerian, and 32 in Japanese SZ patients. These 88 de novo mutations exhibited an enrichment in genes encoding proteins related to GO:0051015~actin filament binding (p = 0.0011) using David, and enrichments in GO: 0003774~transport (p = 0.019) and GO:0003729~mRNA binding (p = 0.010) using Amigo. One of these de novo variant was found in CORO1C coding sequence. We studied Coro1c haploinsufficiency in a Coro1c+/- mouse and found defects in the corpus callosum. These results could motivate future studies of the mechanisms surrounding genes encoding proteins involved in transport and the cytoskeleton, with the goal of developing therapeutic intervention strategies for a subset of SZ cases.

Assessment of Mendelian and risk-factor genes in Alzheimer disease: A prospective nationwide clinical utility study and recommendations for genetic screening.

PURPOSE: To assess the likely pathogenic/pathogenic (LP/P) variants rates in Mendelian dementia genes and the moderate-to-strong risk factors rates in patients with Alzheimer disease (AD). METHODS: We included 700 patients in a prospective study and performed exome sequencing. A panel of 28 Mendelian and 6 risk-factor genes was interpreted and returned to patients. We built a framework for risk variant interpretation and risk gradation and assessed the detection rates among early-onset AD (EOAD, age of onset (AOO) ≤65 years, n = 608) depending on AOO and pedigree structure and late-onset AD (66 < AOO < 75, n = 92). RESULTS: Twenty-one patients carried a LP/P variant in a Mendelian gene (all with EOAD, 3.4%), 20 of 21 affected APP, PSEN1, or PSEN2. LP/P variant detection rates in EOAD ranged from 1.7% to 11.6% based on AOO and pedigree structure. Risk factors were found in 69.5% of the remaining 679 patients, including 83 (12.2%) being heterozygotes for rare risk variants, in decreasing order of frequency, in TREM2, ABCA7, ATP8B4, SORL1, and ABCA1, including 5 heterozygotes for multiple rare risk variants, suggesting non-monogenic inheritance, even in some autosomal-dominant-like pedigrees. CONCLUSION: We suggest that genetic screening should be proposed to all EOAD patients and should no longer be prioritized based on pedigree structure.

Multi-ancestry polygenic risk scores for venous thromboembolism.

Venous thromboembolism (VTE) is a significant contributor to morbidity and mortality, with large disparities in incidence rates between Black and White Americans. Polygenic risk scores (PRSs) limited to variants discovered in genome-wide association studies in European-ancestry samples can identify European-ancestry individuals at high risk of VTE. However, there is limited evidence on whether high-dimensional PRS constructed using more sophisticated methods and more diverse training data can enhance the predictive ability and their utility across diverse populations. We developed PRSs for VTE using summary statistics from the International Network against Venous Thrombosis (INVENT) consortium GWAS meta-analyses of European- (71,771 cases and 1,059,740 controls) and African-ancestry samples (7,482 cases and 129,975 controls). We used LDpred2 and PRSCSx to construct ancestry-specific and multi-ancestry PRSs and evaluated their performance in an independent European- (6,261 cases and 88,238 controls) and African-ancestry sample (1,385 cases and 12,569 controls). Multi-ancestry PRSs with weights tuned in European- and African-ancestry samples, respectively, outperformed ancestry-specific PRSs in European- (PRSCSXEUR: AUC=0.61 (0.60, 0.61), PRSCSX_combinedEUR: AUC=0.61 (0.60, 0.62)) and African-ancestry test samples (PRSCSXAFR: AUC=0.58 (0.57, 0.6), PRSCSX_combined AFR: AUC=0.59 (0.57, 0.60)). The highest fifth percentile of the best-performing PRS was associated with 1.9-fold and 1.68-fold increased risk for VTE among European- and African-ancestry subjects, respectively, relative to those in the middle stratum. These findings suggest that the multi-ancestry PRS may be used to identify individuals at highest risk for VTE and provide guidance for the most effective treatment strategy across diverse populations.

Episignatures in practice: independent evaluation of published episignatures for the molecular diagnostics of ten neurodevelopmental disorders.

Variants of uncertain significance (VUS) are a significant issue for the molecular diagnosis of rare diseases. The publication of episignatures as effective biomarkers of certain Mendelian neurodevelopmental disorders has raised hopes to help classify VUS. However, prediction abilities of most published episignatures have not been independently investigated yet, which is a prerequisite for an informed and rigorous use in a diagnostic setting. We generated DNA methylation data from 101 carriers of (likely) pathogenic variants in ten different genes, 57 VUS carriers, and 25 healthy controls. Combining published episignature information and new validation data with a k-nearest-neighbour classifier within a leave-one-out scheme, we provide unbiased specificity and sensitivity estimates for each of the signatures. Our procedure reached 100% specificity, but the sensitivities unexpectedly spanned a very large spectrum. While ATRX, DNMT3A, KMT2D, and NSD1 signatures displayed a 100% sensitivity, CREBBP-RSTS and one of the CHD8 signatures reached <40% sensitivity on our dataset. Remaining Cornelia de Lange syndrome, KMT2A, KDM5C and CHD7 signatures reached 70-100% sensitivity at best with unstable performances, suffering from heterogeneous methylation profiles among cases and rare discordant samples. Our results call for cautiousness and demonstrate that episignatures do not perform equally well. Some signatures are ready for confident use in a diagnostic setting. Yet, it is imperative to characterise the actual validity perimeter and interpretation of each episignature with the help of larger validation sample sizes and in a broader set of episignatures.

Inverse Modulation of Aurora Kinase A and Topoisomerase IIα in Normal and Tumor Breast Cells upon Knockdown of Mitochondrial ASncmtRNA.

Breast cancer is currently the most diagnosed form of cancer and the leading cause of death by cancer among females worldwide. We described the family of long non-coding mitochondrial RNAs (ncmtRNAs), comprised of sense (SncmtRNA) and antisense (ASncmtRNA) members. Knockdown of ASncmtRNAs using antisense oligonucleotides (ASOs) induces proliferative arrest and apoptotic death of tumor cells, but not normal cells, from various tissue origins. In order to study the mechanisms underlying this selectivity, in this study we performed RNAseq in MDA-MB-231 breast cancer cells transfected with ASncmtRNA-specific ASO or control-ASO, or left untransfected. Bioinformatic analysis yielded several differentially expressed cell-cycle-related genes, from which we selected Aurora kinase A (AURKA) and topoisomerase IIα (TOP2A) for RT-qPCR and western blot validation in MDA-MB-231 and MCF7 breast cancer cells, as well as normal breast epithelial cells (HMEC). We observed no clear differences regarding mRNA levels but both proteins were downregulated in tumor cells and upregulated in normal cells. Since these proteins play a role in genomic integrity, this inverse effect of ASncmtRNA knockdown could account for tumor cell downfall whilst protecting normal cells, suggesting this approach could be used for genomic protection under cancer treatment regimens or other scenarios.

A Case-Only Genome-Wide Interaction Study of Smoking and Bladder Cancer Risk: Results from the COBLAnCE Cohort.

Bladder cancer (BC) is the 6th most common cancer worldwide, with tobacco smoking considered as its main risk factor. Accumulating evidence has found associations between genetic variants and the risk of BC. Candidate gene-environment interaction studies have suggested interactions between cigarette smoking and NAT2/GSTM1 gene variants. Our objective was to perform a genome-wide association case-only study using the French national prospective COBLAnCE cohort (COhort to study BLAdder CancEr), focusing on smoking behavior. The COBLAnCE cohort comprises 1800 BC patients enrolled between 2012 and 2018. Peripheral blood samples collected at enrolment were genotyped using the Illumina Global Screening Array with a Multi-Disease drop-in panel. Genotyping data (9,719,614 single nucleotide polymorphisms (SNP)) of 1674, 1283, and 1342 patients were analyzed for smoking status, average tobacco consumption, and age at smoking initiation, respectively. A genome-wide association study (GWAS) was conducted adjusting for gender, age, and genetic principal components. The results suggest new candidate loci (4q22.1, 12p13.1, 16p13.3) interacting with smoking behavior for the risk of BC. Our results need to be validated in other case-control or cohort studies.

Multiancestry analysis of the HLA locus in Alzheimer's and Parkinson's diseases uncovers a shared adaptive immune response mediated by HLA-DRB1*04 subtypes.

Across multiancestry groups, we analyzed Human Leukocyte Antigen (HLA) associations in over 176,000 individuals with Parkinson's disease (PD) and Alzheimer's disease (AD) versus controls. We demonstrate that the two diseases share the same protective association at the HLA locus. HLA-specific fine-mapping showed that hierarchical protective effects of HLA-DRB1*04 subtypes best accounted for the association, strongest with HLA-DRB1*04:04 and HLA-DRB1*04:07, and intermediary with HLA-DRB1*04:01 and HLA-DRB1*04:03. The same signal was associated with decreased neurofibrillary tangles in postmortem brains and was associated with reduced tau levels in cerebrospinal fluid and to a lower extent with increased Aß42. Protective HLA-DRB1*04 subtypes strongly bound the aggregation-prone tau PHF6 sequence, however only when acetylated at a lysine (K311), a common posttranslational modification central to tau aggregation. An HLA-DRB1*04-mediated adaptive immune response decreases PD and AD risks, potentially by acting against tau, offering the possibility of therapeutic avenues.

Polygenic risk of major depressive disorder as a risk factor for venous thromboembolism.

Major depressive disorder (MDD), bipolar disorder (BD), and schizophrenia (SCZ) are associated with an increased risk of cardiovascular diseases, including venous thromboembolism (VTE). The reasons for this are complex and include obesity, smoking, and use of hormones and psychotropic medications. Genetic studies have increasingly provided evidence of the shared genetic risk of psychiatric and cardiometabolic illnesses. This study aimed to determine whether a genetic predisposition to MDD, BD, or SCZ is associated with an increased risk of VTE. Genetic correlations using the largest genome-wide genetic meta-analyses summary statistics for MDD, BD, and SCZ (Psychiatric Genetics Consortium) and a recent genome-wide genetic meta-analysis of VTE (INVENT Consortium) demonstrated a positive association between VTE and MDD but not BD or SCZ. The same summary statistics were used to construct polygenic risk scores for MDD, BD, and SCZ in UK Biobank participants of self-reported White British ancestry. These were assessed for impact on self-reported VTE risk (10 786 cases, 285 124 controls), using logistic regression, in sex-specific and sex-combined analyses. We identified significant positive associations between polygenic risk for MDD and the risk of VTE in men, women, and sex-combined analyses, independent of the known risk factors. Secondary analyses demonstrated that this association was not driven by those with lifetime experience of mental illness. Meta-analyses of individual data from 6 additional independent cohorts replicated the sex-combined association. This report provides evidence for shared biological mechanisms leading to MDD and VTE and suggests that, in the absence of genetic data, a family history of MDD might be considered when assessing the risk of VTE.

Elevated plasma complement factor H related 5 protein is associated with venous thromboembolism.

Venous thromboembolism (VTE) is a common, multi-causal disease with potentially serious short- and long-term complications. In clinical practice, there is a need for improved plasma biomarker-based tools for VTE diagnosis and risk prediction. Here we show, using proteomics profiling to screen plasma from patients with suspected acute VTE, and several case-control studies for VTE, how Complement Factor H Related 5 protein (CFHR5), a regulator of the alternative pathway of complement activation, is a VTE-associated plasma biomarker. In plasma, higher CFHR5 levels are associated with increased thrombin generation potential and recombinant CFHR5 enhanced platelet activation in vitro. GWAS analysis of ~52,000 participants identifies six loci associated with CFHR5 plasma levels, but Mendelian randomization do not demonstrate causality between CFHR5 and VTE. Our results indicate an important role for the regulation of the alternative pathway of complement activation in VTE and that CFHR5 represents a potential diagnostic and/or risk predictive plasma biomarker.

Genome-wide association study of a semicontinuous trait: illustration of the impact of the modeling strategy through the study of Neutrophil Extracellular Traps levels.

Over the last years, there has been a considerable expansion of genome-wide association studies (GWAS) for discovering biological pathways underlying pathological conditions or disease biomarkers. These GWAS are often limited to binary or quantitative traits analyzed through linear or logistic models, respectively. In some situations, the distribution of the outcome may require more complex modeling, such as when the outcome exhibits a semicontinuous distribution characterized by an excess of zero values followed by a non-negative and right-skewed distribution. We here investigate three different modeling for semicontinuous data: Tobit, Negative Binomial and Compound Poisson-Gamma. Using both simulated data and a real GWAS on Neutrophil Extracellular Traps (NETs), an emerging biomarker in immuno-thrombosis, we demonstrate that Compound Poisson-Gamma was the most robust model with respect to low allele frequencies and outliers. This model further identified the MIR155HG locus as significantly (P = 1.4 × 10-8) associated with NETs plasma levels in a sample of 657 participants, a locus recently highlighted to be involved in NETs formation in mice. This work highlights the importance of the modeling strategy for GWAS of a semicontinuous outcome and suggests Compound Poisson-Gamma as an elegant but neglected alternative to Negative Binomial for modeling semicontinuous outcome in the context of genomic investigations.

Genetic factors for differentiated thyroid cancer in French Polynesia: new candidate loci.

Background: Populations of French Polynesia (FP), where France performed atmospheric tests between 1966 and 1974, experience a high incidence of differentiated thyroid cancer (DTC). However, up to now, no sufficiently large study of DTC genetic factors in this population has been performed to reach definitive conclusion. This research aimed to analyze the genetic factors of DTC risk among the native FP populations. Methods: We analyzed more than 300 000 single nucleotide polymorphisms (SNPs) genotyped in 283 DTC cases and 418 matched controls born in FP, most being younger than 15 years old at the time of the first nuclear tests. We analyzed the genetic profile of our cohort to identify population subgroups. We then completed a genome-wide analysis study on the whole population. Results: We identified a specific genetic structure in the FP population reflecting admixture from Asian and European populations. We identified three regions associated with increased DTC risk at 6q24.3, 10p12.2, and 17q21.32. The lead SNPs at these loci showed respective p-values of 1.66 × 10-7, 2.39 × 10-7, and 7.19 × 10-7 and corresponding odds ratios of 2.02, 1.89, and 2.37. Conclusion: Our study results suggest a role of the loci 6q24.3, 10p12.2 and 17q21.32 in DTC risk. However, a whole genome sequencing approach would be better suited to characterize these factors than genotyping with microarray chip designed for the Caucasian population. Moreover, the functional impact of these three new loci needs to be further explored and validated.

WGS Revealed Novel BBS5 Pathogenic Variants, Missed by WES, Causing Ciliary Structure and Function Defects.

Bardet-Biedl syndrome (BBS) is an autosomal recessive ciliopathy that affects multiple organs, leading to retinitis pigmentosa, polydactyly, obesity, renal anomalies, cognitive impairment, and hypogonadism. Until now, biallelic pathogenic variants have been identified in at least 24 genes delineating the genetic heterogeneity of BBS. Among those, BBS5 is a minor contributor to the mutation load and is one of the eight subunits forming the BBSome, a protein complex implied in protein trafficking within the cilia. This study reports on a European BBS5 patient with a severe BBS phenotype. Genetic analysis was performed using multiple next-generation sequencing (NGS) tests (targeted exome, TES and whole exome, WES), and biallelic pathogenic variants could only be identified using whole-genome sequencing (WGS), including a previously missed large deletion of the first exons. Despite the absence of family samples, the biallelic status of the variants was confirmed. The BBS5 protein's impact was confirmed on the patient's cells (presence/absence and size of the cilium) and ciliary function (Sonic Hedgehog pathway). This study highlights the importance of WGS and the challenge of reliable structural variant detection in patients' genetic explorations as well as functional tests to assess a variant's pathogenicity.

Excessive self-grooming, gene dysregulation and imbalance between the striosome and matrix compartments in the striatum of Shank3 mutant mice.

Autism is characterized by atypical social communication and stereotyped behaviors. Mutations in the gene encoding the synaptic scaffolding protein SHANK3 are detected in 1-2% of patients with autism and intellectual disability, but the mechanisms underpinning the symptoms remain largely unknown. Here, we characterized the behavior of Shank3 Δ11/Δ11 mice from 3 to 12 months of age. We observed decreased locomotor activity, increased stereotyped self-grooming and modification of socio-sexual interaction compared to wild-type littermates. We then used RNAseq on four brain regions of the same animals to identify differentially expressed genes (DEGs). DEGs were identified mainly in the striatum and were associated with synaptic transmission (e.g., Grm2, Dlgap1), G-protein-signaling pathways (e.g., Gnal, Prkcg1, and Camk2g), as well as excitation/inhibition balance (e.g., Gad2). Downregulated and upregulated genes were enriched in the gene clusters of medium-sized spiny neurons expressing the dopamine 1 (D1-MSN) and the dopamine 2 receptor (D2-MSN), respectively. Several DEGs (Cnr1, Gnal, Gad2, and Drd4) were reported as striosome markers. By studying the distribution of the glutamate decarboxylase GAD65, encoded by Gad2, we showed that the striosome compartment of Shank3 Δ11/Δ11 mice was enlarged and displayed much higher expression of GAD65 compared to wild-type mice. Altogether, these results indicate altered gene expression in the striatum of Shank3-deficient mice and strongly suggest, for the first time, that the excessive self-grooming of these mice is related to an imbalance in the striatal striosome and matrix compartments.

High diagnostic potential of short and long read genome sequencing with transcriptome analysis in exome-negative developmental disorders.

Exome sequencing (ES) has become the method of choice for diagnosing rare diseases, while the availability of short-read genome sequencing (SR-GS) in a medical setting is increasing. In addition, new sequencing technologies, such as long-read genome sequencing (LR-GS) and transcriptome sequencing, are being increasingly used. However, the contribution of these techniques compared to widely used ES is not well established, particularly in regards to the analysis of non-coding regions. In a pilot study of five probands affected by an undiagnosed neurodevelopmental disorder, we performed trio-based short-read GS and long-read GS as well as case-only peripheral blood transcriptome sequencing. We identified three new genetic diagnoses, none of which affected the coding regions. More specifically, LR-GS identified a balanced inversion in NSD1, highlighting a rare mechanism of Sotos syndrome. SR-GS identified a homozygous deep intronic variant of KLHL7 resulting in a neoexon inclusion, and a de novo mosaic intronic 22-bp deletion in KMT2D, leading to the diagnosis of Perching and Kabuki syndromes, respectively. All three variants had a significant effect on the transcriptome, which showed decreased gene expression, mono-allelic expression and splicing defects, respectively, further validating the effect of these variants. Overall, in undiagnosed patients, the combination of short and long read GS allowed the detection of cryptic variations not or barely detectable by ES, making it a highly sensitive method at the cost of more complex bioinformatics approaches. Transcriptome sequencing is a valuable complement for the functional validation of variations, particularly in the non-coding genome.